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We present a careful recasting of first-principles band structure calculations for MgB2 in a non- 
orthogonal sp-tight-binding (TB) basis. Our TB results almost exactly reproduce our full potential 
linearized augmented plane wave results for the energy bands, the densities of states and the total 
energies. Our procedure generates transferable Slater-Koster parameters which should be useful for 
other studies of this important material. 



The recent discovery of superconductivity in MgB2EJ 
has created great interest in the study of this material, 
both to understand the mechanism of superconductiv- 
ity and to explore other properties of MgB2 and related 
materials. Intensive, research has beierjt. carried out both 
by experimentalistsErcl and theorists.oia There have been 
several studies of the electronic structure of MgB2 in- 
cluding total energy, band structure and phonon spectra 
calculations as well as evaluations of the electron-phonon 
coupling, which seems to have emerged as the prime can- 
didate for explaining the superconducting behavior. 

In this paper we present a highly accurate tight- 
binding (TB) description of the band structure and total 
energy of MgB2. While there have been TB interpreta- 
tions of the electronic structure of MgB2 in the literature, 
a realistic recasting of the details of the first-principles 
electronic structure calculations is larking. Our approach 
follows the NRL-TB methodology,E30 which is based 
on deriving a non-orthogonal TB Hamiltonian by fit- 
ting to both the total energy and energy band results 
of a first-principles fuLLpatential Linearized Augmented 
Plane Wave (LAPW)B'B calculation using the Hedin- 
Lundqvist paraipetrization of the Local Density Approx- 
imation (LDA).EJ We first performed detailed LAPW 
calculations for MgB2 in its ground state (AIB2) struc- 
ture, varying c and a, thus determining the LDA equi- 
librium volume. It was necessary to perform 17 indepen- 
dent LAPW calculations over a large range of volumes 
and c/a ratios. Our LAPW equilibrium parameters are 
c = 6.55 a.u. and a = 5.75 a.u., as compared to the ex- 
perimental values of c = 6.66 a.u. and a = 5.83 a.u. As 
is usually the case the LDA underestimates the experi- 
mental values, here by about 1.5%. 

All the above results, i.e., 17 values of the total energy 
and the energy bands for 76 k-points in the irreducible 
hexagonal Brillouin zone, were used as a database to de- 
termine the parameters of our tight-binding Hamiltonian. 
According to the NRL-TB scheme the on-site parame- 
ters depend on the density of the neighboring atoms and 
the hopping integrals have a polynomial dependence that 
extends to at least the third nearest-neighbor distance. 
Our basis included the s and p orbitals on both Mg and B 
in a non-orthogonal two-center representation. A wave- 
function analysis of our LAPW results shows that the 
bands up to the Fermi level, ep, are strongly dominated 
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FIG. 1. The band structure of MgB2 in the AIB2 struc- 
ture at the theoretical equilibrium volume, as determined 
by the full-potential LAPW method (solid lines) and our 
tight-binding parametrization (dashed lines). The Fermi level 
is at e = 0. 

by the B p states with very little contribution from the 
Mg ions. It turns out, however, that an accurate TB fit 
including only the B orbitals is impossible, and therefore 
the Mg s and p orbitals were included in the fit. Fur- 
thermore, to obtain a highly accurate fit it was essential 
to block diagonalize the Hamiltonian at the high sym- 
metry points r. A, L, K and H. We find that at a given 
set of lattice parameters (c,a) we can perfectly reproduce 
the energy bands of MgB2. A comparison is shown in 
Fig. 0, where the solid and broken lines represent the 
LAPW and TB bands, respectively, at the LDA values 
of the equilibrium lattice parameters. The TB bands are 
in excellent agreement with the LAPW bands, including 
the two-dimensional B-cr band in the T —> A direction 
just above ep, which has becji, identified as hole band 
controlling superconductivityElQ The RMS fitting error 
is 2 mRy for the total energy, and close to 10 mRy for 
the first five bands. Beyond the fifth band our fit is not 
as accurate, as the Mg d-bands, which are not included in 
our Hamiltonian, come into play. The values of our TB 
parameters are givep in Table m following the notation 
of Bernstein et al£3 In this table we also show, for the 
convenience of the reader, the actual Slater-Koster pa- 
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FIG. 2. The electronic density of states (DOS) of MgB2 
in the AIB2 structure at the theoretical equilibrium volume, 
comparing the total DOS as determined by the full-potential 
LAPW method (upper solid line) and our tight-binding 
parametrization (upper dashed line), and the partial sin- 
gle-atom B-p decomposition (lower lines). The LAPW result 
decomposition was determined inside the muffin-tin and then 
scaled by a factor of 2.37 (see text). 

rameters for three or four nearest neighbors determined 
from our formulas for the specific LDA equilibrium values 
of the lattice constants. 

In Fig. H we show a comparison of TB and LAPW 
densities of states. There is an excellent agreement in 
both the total DOS and the B p-\ike DOS. To facihtate 
the comparison we have normalized the muffin-tin de- 
composed LAPW values so that the contributions from 
the angular momentum components add up to the total 
DOS, as is the case in the TB. For the Boron states this 
amounted to multiplying the decomposed values by 2.37. 
The B s-components of the DOS have their strongest 
presence at the bottom of the valence band, from -0.8 Ry 
to -0.6 Ry on our scale. They are much smaller than the 
p-like DOS, so we chose not to include them in Fig. ||. 
Additionally, we have omitted the Mg p-like DOS, which 
is also small below ep, although it becomes significant 
above ep- Our TB value of the total DOS at ep is 
N{eF) = 0.69 states/eV, which is almost identical to that 
found from our direct LAPW calculation. This value of 
N^ep) corresponds to the LDA equilibrium volume and 
is slightly smaller thanJhfi value of 0.71 states/eV re- 
ported by other workerajl3 at the experimental jiplume. 
Using our value of N{ep) and the measured valucH of the 
specific heat coefficient 7 we infer a value of A = 0.65 
which is consistent with the high superconducting tran- 
sition temperature in MgB2. It should also be noted that 
the B-p states contribute 81% of the DOS at ep. 

Our TB-Hamiltonian also provides an accurate de- 
scription of the energetics of MgB2, which is expected 
to be very useful for other theoretical studies. We have 
further tested our parameters by computing the TB equi- 
librium structure. We find an equilibrium of c = 6.66 a.u. 



and a = 5.79 a.u., in good agreement with the LAPW re- 
sult. At c/a = 1.14 (the experimental value), we deduce 
a bulk modulus of B = 165 GPa which in good agree- 
ment with the experimental value of 120GPa and with 
the. calculated value of 147 GPa reported by Bohnen et 
alE 

The TB parameters presented in this paper give a very 
accurate description of the band structure of MgB2. The 
availability of this Hamiltonian should motivate the cal- 
culation of other properties of this important material. 

We thank I. I. Mazin and D. J. Singh for useful dis- 
cussions. This work was supported by the U. S. Office of 
Naval Research. The development of the tight-binding 
codes was supported in part by the U. S. Department 
of Defense Common HPC Software Support Initiative 
(CHSSI). 
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TABLE I. Tight- binding parameters for MgB2, generated following the methods of Mehl and Papaconstantopoulostil and 
Bernstein et aLtJ. Also shown are the generated Slater-Koster tight-binding parameters for the nearest neighbors at the LDA 
equilibrium lattice constants, a = 5.75 a.u. and c = 6.53 a.u. On-site energies are generated from the "densities" of like atoms, 
that is, the Mg on-site parameters come from the density of Mg atomfc-iand the B on-site parameters from the density of B 
atoms. F{R) is the cutoff function from equation (2) of Berstein et alxB, with Rc = 12.5 a.u. and Lc — 0.5a.u. All energies 
are in Rydbergs, all distances in a.u. 

Mg-Mg Interactions 



On-site Parameters (A = 0.93961 a.ur'^^'^) 



e 


Q-e 


Pi 


7« 


Xe 


LDA Equilibrium Values 


s 


0.02169 


-0.25368 


-0.04017 


19.84215 






0.03516 




p 


0.39868 


-0.22303 


1.35834 


53.36624 






0.52322 








Hopping Terms: Hui^{R) 


= + bii 


pi? + Cui^R'') 


eM-9u'^R)F{R) 










Cet'ij. 






a 


c 


\/3a 


H(ssa) 


5715.097 


-310.8836 


-182.0526 


1.35579 




-0.05372 


-0.02495 


-0.00009 


H{spa) 


5704288. 


541286.7 


-387450.5 


1.84506 




-0.01259 


-0.00161 


-0.00000 


H{ppa) 


-1920.935 


498.3775 


-22.58129 


1.12482 




0.13720 


0.09557 


0.00141 


HijppTi) 


2000.513 


-739.8181 


70.26517 


1.13170 




0.04414 


0.03861 


0.00241 






Overlap Terms: Sui i^i_{R) = {5t 


e' + tu'tj.R + 


lU'tJ.R'' + 


i?^)exp(-u^,,^i?)f(7?) 




Su' fj. 




<]ee'fj. 








a 


c 


\/3a 


S{ssa) 


1.04886 


-1.27181 


0.55382 


1.01629 




0.18512 


0.12683 


0.00767 


S{spa) 


0.41781 


0.03630 


-0.00873 


0.63396 




0.19274 


0.13373 


-0.00815 


S{ppa) 


-24.36368 


0.17541 


0.41661 


1.07340 




-0.07174 


-0.01867 


0.00101 


S{pptt) 


-68.95974 


5.97517 


2.55826 


1.19338 




0.08007 


0.04731 


0.00088 



B-B Interactions 



On-site Parameters (A = 0.79205 a.M."^^^) 



£ 


at 


111 


Je 


xe 


LDA Equilibrium Values 


s 


-0.16521 


-0.00022 


0.02579 


0.09088 




-0.09356 






P 


0.38802 


0.00060 


0.00566 


0.01918 




0.40383 










Hopping Terms 


'■ Hiiifj_{R) 


— + bu'^R + cuifj^R''] 


^M-gu'uR)nR) 










bu'fi 


Cee'ij. 


gw/j. 


a/^/3 


a 


c 


2a/x/3 


H{ssa) 


-7.31550 


2.09241 


-0.23379 


0.85573 


-0.25908 


-0.04471 


-0.03032 


-0.02881 


H{spa) 


-146.6733 


64.75572 


-8.57386 


1.21983 


-0.18743 


-0.01112 


-0.00538 


- 0.00484 


Hippo-) 


-296.2214 


128.0942 


10.42990 


1.17220 


0.14703 


0.03537 


0.01210 


0.01030 


H(ppTi) 


167.1287 


-84.95581 


9.40729 


1.16057 


-0.12834 


-0.00448 


0.00204 


0.00232 






Overlap Terms: Set 


'm(^) = iS> 


e> + tei'i^R + Qee'fiR'' + f'li'pi 


i?^)exp(-<,^i?)F(i?) 






Su'fj, 




<lU' fj. 






a/^/3 


a 


c 


2a/x/3 


S{ssa) 


0.08974 


-0.05865 


0.00446 


0.60130 


0.24535 


0.05320 


0.03092 


0.02870 


S{spa) 


14.02893 


-2.43293 


0.33914 


1.27302 


0.14823 


0.00581 


0.00208 


0.00181 


S{ppa) 


-60.70629 


-0.47590 


3.29498 


1.32165 


-0.25833 


0.01141 


0.00558 


0.00497 


S{ppTv) 


18.98764 


6.13369 


-3.44742 


1.41889 


0.00689 


-0.00321 


-0.00112 


-0.00096 


Mg-B Interactions 






Hopping Terms 


: Hu'tj.{R) 


= (a«v + bu'^^R + Cu'^^R^) exp{-gfi,,^^R)F{R) 










bte'fx 


Cte'ij. 








a? -f ic2 




H(ssa) 


-15.40626 


8.92332 


-2.25890 


1.06263 


-0.11887 


-0.01709 


-0.00257 




H{spa) 


-22.65145 


5.35089 


-0.60679 


1.03205 


-0.07642 


-0.00612 


-0.00093 




H{ppa) 


98.38228 


-45.01479 


6.05711 


1.20823 


0.02245 


0.00196 


0.00019 




H{pptt) 


-94.47230 


33.60639 


-4.25418 


1.21106 


-0.03269 


-0.00152 


-0.00013 




H{psa) 


7.80580 


1.71300 


-0.22442 


1.03201 


0.07662 


0.00308 


0.00015 








Overlap Terms 


'■ Su'ti{R) 


= {tu't^ + q 


u'fj.R + rw^R'') 


eM-<i',R)F{R) 






Su'fj, 










Via2 + lc2 




a? ^ ic2 




S{ssa) 


1.74820 


0.13546 


0.07434 


0.82425 


0.16873 


0.04450 


0.01278 




S{apa) 


15.27243 


3.72217 


0.27051 


1.09173 


0.14959 


0.00848 


0.00081 




S{ppa) 


-4.51769 


-4.96427 


0.71842 


0.95339 


-0.17503 


-0.00229 


0.00188 




S{ppTv) 


846.58108 


265.43883 


10.03535 


1.30600 


0.06108 


0.00187 


0.00007 




S{psa) 


-2.81156 


0.20700 


-0.40404 


0.92459 


-0.19812 


-0.04168 


-0.00940 
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